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Abstract 

A new finite atlas of overlapping balanced canonical forms for multivariate discrete-time loss- 
less systems is presented. The canonical forms have the property that the controllability matrix is 
positive upper triangular up to a suitable permutation of its columns. This is a generalization of a 
similar balanced canonical form for continuous-time lossless systems. It is shown that this atlas is 
in fact a finite sub-atlas of the infinite atlas of overlapping balanced canonical forms for lossless 
systems that is associated with the tangential Schur algorithm; such canonical forms satisfy certain 
interpolation conditions on a corresponding sequence of lossless transfer matrices. The connection 
between these balanced canonical forms for lossless systems and the tangential Schur algorithm 
for lossless systems is a generalization of the same connection in the SISO case that was noted 
before. The results are directly applicable to obtain a finite sub-atlas of multivariate input-normal 
canonical forms for stable linear systems of given fixed order, which is minimal in the sense that 
no chart can be left out of the atlas without losing the property that the atlas covers the manifold. 

Keywords: Lossless systems, input normal forms, output normal forms, balanced canonical forms, 
model reduction, MIMO systems, tangential Schur algorithm. 



1 Introduction 

In linear systems theory there has been a longstanding program in developing balanced realizations, 
balanced canonical forms and associated parameterizations for stable linear systems and for various 
other classes of linear systems. The classical Gramian based concept of balancing as introduced by 
Moore, see 1 10], applies to stable systems and allows one to develop parameterizations in which system 
stability is a built-in property. One of the motivations for the interest in balancing is that it leads to a 
simple method for model order reduction, namely by truncation of (the last entries of) the state vector. 

However truncation does not always lead to a minimal system. Therefore there has been research 
into balanced canonical forms which do have the property that truncation of the last entries in the state 
vector leads to a minimal system. For continuous-time systems this has led to the original balanced 
canonical form of Ober (see [11]) and to the new balanced canonical form of Hanzon (see lO; see 
also [12]). This last balanced canonical form is based on the idea that if the controllability matrix is 
positive upper triangular (i.e., the controllability matrix forms an upper triangular matrix with positive 
entries on the pivot positions), then truncation of the last entries of the state vector leads again to a 



system with a positive upper triangular controllability matrix, hence is controllable. Because this is in 
the balanced continuous-time case, the controllability property here implies that the resulting system 
is again minimal and balanced. 

To use similar ideas to build overlapping balanced canonical forms is more involved. For continuous- 
time lossless systems, which form the key to these problems, a generalization of positive upper trian- 
gular matrices is used in {€]. The idea used there is that it suffices if a column permutation of the 
controllability matrix is positive upper triangular. Under certain circumstances there will exist an as- 
sociated column permutation (we also speak of a shuffle of columns in this context) of the so-called 
realization matrix, which allows one to proceed with the construction. 

In the case of discrete-time systems the situation is somewhat more complicated because it is 
known that starting from a balanced realization, truncation of the state vector will normally not lead to 
a balanced state-space system. In the case of SISO lossless discrete-time systems a balanced canonical 
form with a simple positive upper triangular controllability matrix was presented in jSl. Also the 
possibilities for model reduction by truncation, combined with a correction of some sort to arrive at a 
balanced realization of a lossless system, are discussed there. 

In the current paper we treat the case of MIMO lossless discrete-time systems. We present overlap- 
ping balanced canonical forms which have the property that the corresponding controllability matrix 
is positive upper triangular, up to a column permutation. In this sense it is close to the corresponding 
results in [6]; however, here a generalization is presented which simplifies the presentation and which 
can, as a spin-off, also be used in the continuous-time case. The precise relation with the approach 
taken in 16] will be made clear. The results on the relation between a specific triangular pivot structure 
in controllable pairs, which we call "staircase forms", and an associated triangular pivot structure in 
the controllability matrices are also of interest outside the context of lossless systems. 

In ISl a connection was shown between the balanced canonical forms there presented and the Schur 
algorithm for scalar lossless discrete-time transfer functions. In Q it is shown how the parameteriza- 
tions for multivariable rational lossless transfer matrices by Schur parameters, based on the so-called 
tangential Schur algorithm, can likewise be lifted into parameterizations by Schur parameters of bal- 
anced state-space canonical forms of lossless systems. One of the main results of the current paper is 
to show how the atlas of overlapping balanced canonical forms presented in this paper can be obtained 
as a finite sub-atlas of the infinite atlas of overlapping balanced canonical forms corresponding to the 
tangential Schur algorithm. In fact, a certain well-specified choice of so-called direction vectors in the 
tangential Schur algorithm leads to the balanced canonical forms presented here. 

Although a generalization of the results of this paper to the case of complex-valued systems is 
straightforward, we shall restrict the discussion to the case of real-valued systems only for ease of 
presentation. 

2 Preliminaries 

2.1 State space systems and realization theory 

Consider a hnear time-invariant state-space system in discrete time with m inputs and m outputs: 



with t G 7^, xt E M" for some nonnegative integer n (the state space dimension), ut G M*" and 
yt G M™. The matrices A, B, C and D with real-valued entries are of compatible sizes: n x n, 



xt+i = Axt + But 
yt = Cxt + Dut, 



(1) 
(2) 
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n X m, m X n and m x m, respectively. The corresponding transfer matrix of this system is given by 
G{z) = D + C{zln — A)~^B, which is an m x m matrix with rational functions as its entries. The 
controllability matrix K and the observability matrix O associated with this system are defined as the 
block-partitioned matrices 



K= [B,AB,...,A^-'B] 



O 



C 

CA 



CA 



n-l 



(3) 



The system (or its input pair {A, B)) is called controllable if K has full row rank n and the system 
(or its output pair (C, A)) is called observable if O has full column rank n. Minimality holds iff both 
controllability and observability hold, which holds iff the McMillan degree of G{z) is equal to n. 

To any such state-space system we associate the following (square) block-partitioned matrix R, 
which we call the realization matrix: 

DC' 
B A 



R 



(4) 



The matrix R, its nx {m + n) sub-matrix [-B, and the associated n x nm controllability matrix K 
will all play an important role in the sequel. 



2.2 Stability and balancing 

Let {A, B,C, D) be some state space realization of a transfer matrix G{z). If the eigenvalues of 
A all belong to the open unit disk in the complex plane, then the matrix A is called (discrete-time) 
asymptotically stable, and [A, B,C, D) is an asymptotically stable realization of G{z). (For more 
details on state-space realization theory, see e.g. [9|.) 

If {A, B, C, D) is an asymptotically stable realization, then the controllability Gramian Wc and the 
observability Gramian Wq are well defined as the exponentially convergent series 

oo 

Wc = ^A''BB^{A^)'', (5) 

k=0 

oo 

Wo = ^{A'^)''G^CA''. (6) 

k=0 

These Gramians are characterized as the unique (and positive semi-definite) solutions of the respective 
Lyapunov-Stein equations 

Wc-AWcA'^ = BB^, (7) 
Wo-A^WoA = G^G. (8) 

A minimal and asymptotically stable state-space realization (^4, B, G, D) of a transfer matrix is called 
balanced if its controllability and observability Gramians Wc and Wo are both diagonal and equal. 
Minimality implies that Wc and Wo are non-singular, hence positive definite. Any minimal and asymp- 
totically stable realization {A, B,G, D) is similar to a balanced realization, meaning that there exists a 
nonsingular state space transformation matrix T which makes the realization {TAT^^,TB, GT^^,D) 
into a balanced realization. 



3 



A system is called input-normal if Wc = In and it is called output-normal if Wo = In- Balanced 
realizations are directly related to input-normal and output-normal realizations, respectively, by diago- 
nal state space transformations. The property of input-normality (resp. output-normality) is preserved 
under orthogonal state space transformations. 



2.3 Lossless systems, balanced realizations and the tangential Schur algorithm 



A discrete-time system is called lossless if it is stable and its m x m transfer matrix G{z) is unitary for 
all complex z with \z\ = 1. It is well-known (cf., e.g.. Proposition 3.2 in iH and the references given 

D C 

is a balanced realization matrix of a lossless system if and only if R is an 



there) that R 



B A 



orthogonal matrix and A is asymptotically stable. It then holds that Wc = Wo = In- For a further 
background on lossless systems, see e.g. f3l. 

In [7] an atlas of overlapping balanced canonical forms for lossless discrete-time systems of order 
n is presented. Also, a closely related atlas is given for (controllable) input-normal pairs {A, B) by 
considering the quotient space with respect to the orthogonal group. Each of these balanced canonical 
forms is then characterized (in the real case) by a fixed sequence of n interpolation points Wk € M, 
\wk\ < 1, = 1, . . . , n, and a fixed sequence of n normalized direction vectors Uk € M"^, \\uk\\ = 1, 
k = 1, . . . , n (which are not to be confused with the input signal applied to a system). Here we shall 
consider the case Wk = 0, k = 1, . . . ,n, hence each balanced canonical form that we consider is 
determined entirely by the choice of direction vectors. Each such balanced canonical form for input- 
normal pairs {A, B) is then parameterized by a sequence of n Schur vectors Vk € M™, with \\vk\\ < 1 
for all A; = 1, . . . ,n. For lossless systems the parameterization also involves an additional m x m 
orthogonal matrix Dq. 

In fact, the realization matrix R in this set-up can be written as an orthogonal matrix product: 



R — ^n^n-i • ■ ■ Ti-RoAf • • • A^, 



(9) 



where for = 1, . . . , n: 





In~k 







Vk 










Ik-l 




" In-k 





Afc = 


Uk 










Ik- 



with an (m + 1) x (m + 1) orthogonal matrix block Vk given by 

Vk 



Vk 



\/l - ll-Vfcl 



(for Vk = Oit holds that Vk 



I™ 

1 

and an (m + 1) x (m + 1) orthogonal matrix block Uk given by 



which makes that Vk depends smoothly on the entries of Vk) 



Uk 



Uk 





- Ukul 

,T 
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and furthermore an (n + m) x (n + m) orthogonal matrix Rq given by 



in which Dq is m x m orthogonal. 

The interpolation conditions attain the form Gk{w^^)uk = v^, where Gk{z) denotes the transfer 
function associated with the k-th order lossless system for which the right lower (m + /c) x (m + k) 
sub-matiix of Rf^ = Tf:Tk_i ■ ■ ■ TiRqAJ ■ ■ ■ A^^A^ is a realization matrix. In the present situation 
with Wk = it follows that G^iw^^) = Gk{oo) = D^., so that the interpolation conditions can be 
written as 

DkUk = Vk 

where {A^, Bk,Ck, D^) denotes the corresponding state-space realization of the fc-th order lossless 
function Gk{z). 

Note that here we consider the real case with real direction vectors and real Schur parameter 
vectors. Note further that Rq, Fi, . . . , r„ and Ai, . . . , A„ are all orthogonal matrices. It is important to 
note that the orthogonal matrix product r„r„_i • • • TiRq in fact forms a positive m-upper Hessenberg 
matrix, i.e. an (m+n) x (m+n) matrix of which the m-th sub-diagonal has positive entries only and of 
which the last n — 1 sub-diagonals are all zero. It also follows that if the direction vectors ui, . . . , n„ 
are taken to be standard basis vectors, then the matrix product Af A|^ • • • A^ yields a permutation 
matrix. Hence in that case the balanced realization matrix R is obtained as a column permutation of an 
orthogonal positive m-upper Hessenberg matrix. 



3 Triangular structures in controllable pairs and their controllability 
matrices 

It is not difficult to see that if the realization matrix R is positive m-upper Hessenberg, then (i) the first 
n columns of the partitioned nx {m + n) matrix [B, A] form a positive upper triangular matrix, i.e. an 
upper triangular matrix with only positive entries on the main diagonal, and (ii) the first n columns of 
the corresponding controllability matrix K = [B, AB, . . . , A^^^B] also form a positive upper triangu- 
lar matrix. (A matrix with this property is called a simple positive upper triangular matrix.) Therefore 
the realization is controllable. In the discrete-time lossless case, if R is orthogonal, controllability 
implies that A is asymptotically stable which in turn implies that the realization is minimal. 

A balanced realization of a lossless system is determined up to an arbitrary orthogonal change of 
basis of the state space. The effect of such a change of basis on the controllability matrix is that it is 
pre-multiplied with an orthogonal matrix. Now it is well-known that any nonsingular square matrix 
can be written as the product of an orthogonal matrix and a positive upper triangular matrix in a unique 
way (in numerical linear algebra this is known as the QR-decomposition). If the first n columns of 
the controllability matrix are linearly independent then a unique orthogonal state-space isomorphism 
exists which transforms the first n columns of the controllability matrix into a positive upper triangular 
matrix. This determines a unique local balanced canonical form for lossless systems. In the SISO case 
it is in fact a global balanced canonical form and it is presented and investigated in [8|. 

In the MIMO case, the canonical form does not apply to systems which have a non-generic Kro- 
necker structure. This is why this is a local canonical form. In order to see how the concept of requiring 
the first n columns of the controllability matrix K to be positive upper triangular can be generahzed 



5 



to obtain an atlas of local canonical forms in the MIMO case, we will consider the relation between 
triangular structures in the partitioned matrix [B, A] and triangular structures in the corresponding 
controUabiUty matrix K = [B, AB, . . . , A^'^B]. The following definitions will turn out to be useful. 

Definition 3.1 Let n be a fixed positive integer. Consider a vector v G M". 

(a) The vector v is called a pivot vector with a pivot at position k, or a pivot-k vector for short, if 
G {1, . . . , n} is an integer for which the entry v(k) is strictly positive and the entries v{j) with j > k 
are all zero. 

(b) The vector v is called a positive pivot-fc vector if it is a pivot-k vector for which in addition the 
entries v{j) with j < k are all strictly positive too. 

Definition 3.2 For given positive integers n and r, consider a mapping J : {1, . . . ,n} —> {0,1, . . . ,r} 
which is written in short-hand notation as J = {ji,j2, ■ ■ ■ 

(a) Associated with J, the mapping : Dj — t- TZj is defined as the restriction of J to Dj" which is 
the largest subset of {1, ... ,n} on which J is nonzero; the co-domain TZj is the corresponding range 
of positive values occurring as images under J. 

(b) The mapping J is called a pivot structure if is a bijection. Then the inverse of is denoted 
by and the extended mapping Q : {1, . . . ,r} — > {0, 1, . . . , n} is written in short-hand notation as 
Q = {q1jQ2, ■ ■ ■ ,Qr} and defined by: qk = Q~^{k)for k G TZj, and Qk = otherwise. 

(c) An nx r matrix M is said to have a pivot structure J if for each k G Vj it holds that column jk of 
M is a pivot-k vector (Equivalently, each column i ofM is a pivot-qi vector, where 'a pivot-0 vector' 
is synonymous to 'not a pivot vector'.) 

(d) A pivot structure J is called a full pivot structure ifDj = {1, . . . , n}. 

Example. Let n = 5 and r = 8. Consider the mapping J : {1, . . . , 5} -> {0, 1, . . . , 8} given by 
J = {ii ) J2, J3, j4, is} = {7, 1, 5, 3, 6}. It follows that the domain and co-domain of J+ are given by 
Vj = {1, 2, 3, 4, 5} and TZj = {1, 3, 5, 6, 7}, respectively. Note that J"*" is a bijection, so that J is 
a pivot structure. Since Vj = {1, . . . , n} it holds that J defines a. full pivot structure. The mapping 
Q, which extends the inverse mapping of J"*", is given by: Q = {qi, q2, 13, Qa, 15, qe, qv : Qs} = 
{2, 0, 4, 0, 3, 5, 1, 0}. Any 5x8 matrix M which has the full pivot structure J is of the following form: 



* 


* 


* 


* 


* 


* 


+ 


* 


+ 


* 


* 


* 


* 


* 





* 





* 


* 


* 


+ 


* 





* 





* 


+ 


* 





* 





* 





* 





* 





+ 





* 



where the entries denoted by * are allowed to have an arbitrary value and the entries denoted by + are 
required to be (strictly) positive. Note that J addresses the entries denoted by + for each row and Q 
specifies the same entries for each column. 

The construction of Q from a given pivot structure J induces a mapping : J ^ Q. From the fact 
that J+ and are each others inverse, is not difficult to see that Tj-^n provides the inverse of r„ 
The sets {{k,jk) \ k G T)^} and {{it, I) \ I G "T^j} obviously coincide: the mappings J and Q both 
serve to specify the same set of pivot locations for an n x r matrix, satisfying the rule that in each row 
and in each column of that matrix at most one pivot location is selected. The mapping J specifies these 
pivot locations in a row-oriented fashion, the mapping Q in a column-oriented fashion. 
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For a full pivot structure it holds that n < r. If J = {ji,j2, ■ ■ ■ , jn} is a full pivot structure for 
an n X r matrix M, then the ordered selection of columns ji,j2,---- jn from M constitutes a positive 
upper triangular nxn sub-matrix. In this way, positive upper triangularity is generalized by the concept 
of a full pivot structure. 

As explained before, if a block-partitioned n x (m + n) matrix [B, A] is simple positive upper 
triangular (i.e., it has the full pivot structure J with = fc for = 1, . . . , n) then the associated (finite 
or infinite) controllability matrix K = [B, AB, A'^B,...] also is simple positive upper triangular. 
We now proceed to investigate the question which full pivot structures for [B, A] induce fuU pivot 
structures for K. Conversely, it is of interest to determine which full pivot structures for K are induced 
by full pivot structures for [B, The latter question is more involved and it is studied in detail in the 
following section. Here we address the former question for which the following definition is useful. 

Definition 3.3 Let m and n be given positive integers. 

(a) A pivot structure F for an n x n matrix A is called a staircase form /or A if is monotonically 
increasing having the range TZ'^ = {1, 2, . . . Here pA denotes the number of pivots, i.e. the 
number of elements in Dp. 

(b) A pivot structure J = {ji, . . . ,jn}for an n x [m + n) block-partitioned matrix [B, A\ induces a 
pivot structure P = {pi , . . . , Pn}for the matrix A as given by p^ = m.ax{jk — m, 0}for k = 1, . . . ,n. 

(c) A full pivot structure J = {ji ,j2,..., jn\ for an n x (m + n) block-partitioned matrix [B, A] is 
called an admissible pivot structure /or [B, A] if it holds that: (i) B has a pivot-1 vector, i.e. I < ji < 
m, and (ii) the pivot structure P induced by J constitutes a staircase form for A. 

Of course, a pivot structure J = {ji, . . . , j„} for an n x [m + n) block-partitioned matrix [B,^] 
also induces a pivot structure for the matrix B. For several purposes, the induced pivot structures 
for A and B are more conveniently described in terms of the associated column-oriented description 
Q = {QI: ■ ■ ■ ) Qm+n} for [^j For the matrix A it holds that the associated column-oriented pivot 
structure S = {si, . . . , s„} satisfies = g„t+fc for all /c = 1, . . . , n. For the matrix B the associated 
column-oriented pivot structure is the restriction of Q to the domain {!,..., m}, simply described by 
the sequence {gi, . . . , qm}- 

Example. Let m = 4, n = 6 and consider the full pivot structure J = {3,1,5,6,4,7} for the 
n X [m + n) partitioned matrix [i?,v4]. The corresponding column-oriented description is given by 
Q = {2, 0, 1, 5, 3, 4, 6, 0, 0, 0}. The matrix [B, A] therefore has the form: 

-|- * * * * * 
+ * * * * 
* * * * 
-|- * * * 

The induced pivot structure for the matrix A is given by P = {pi,P2,P^-,Pi,Pb:P&} = {0,0,1,2,0,3}, 
which follows from p^ = max{jfc — 4, 0} for A; = 1,2,..., 6. The associated column-oriented de- 
scription is then given by S = {si, S2, S3, ^4, S5, se} = {3, 4, 6, 0, 0, 0} = {^5, ge, 97, 98, 99, gio}- The 
function P"*" is given by the pairs (3,1), (4,2) and (6,3); the inverse is given by (1,3), (2,4) 
and (3, 6). Clearly, P+ is monotonically strictly increasing (and equivalently is monotonically 
strictly increasing) so that P is a staircase form for A. This is clearly illustrated by the pattern con- 
stituted by the entries denoted by -I- in the matrix A above. Also, the matrix B has a pivot-1 column 



[B,A\ 



+ 







+ 









* 
* 
* 

* 

+ 
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as its third column. Therefore, J constitutes an admissible pivot structure for [-B,^]. Note that the 
column-oriented description of the pivot structure for the matrix B follows from the restriction of Q 

as: {gi, 92, 93, 94} = {2,0, 1,5}. 

Note that an admissible pivot structure J for [B, A] is totally determined by the induced pivot struc- 
ture for B. In that case, the pivot structure S for the matrix A having a staircase form is given by 
{si, . . . , s„} = {qm+i, • • • , gm+PA' 0' • • • , 0} where the subsequence {g^+i, • • • , Qtu+pa} is positive 
and monotonically increasing, consisting of the elements of {1,2,..., n} not occurring in {qi, . . . , Qm}- 
For admissibility, the only condition on the column-oriented pivot structure {qi, . . . , q^} for B is that 
1 occurs in this sequence. 

If v is a pivot-A; vector and J is admissible, then the staircase structure of A implies that w = Av is a 
pivot-Sfe vector. For this reason, the function S will be called the successor function. (For convenience 
we also define 5(0) = and we recall that the terminology 'a pivot-0 vector' is synonymous to 'not 
a pivot vector'.) The sequence of pivot positions for the vectors Av, A?v, A^v, ... is then given by 
k, S{k), S'^{k), S'^{k), . . .. Conversely, the induced pivot structure P for A is called the predecessor 
function (here we also introduce P{0) = 0). Recall that 5+ and P+ are each others inverse. 

We have the following result. 

Theorem 3.4 Let m and n be given positive integers. 

(a) If J is an admissible pivot structure for an n x [m + n) block-partitioned matrix [B,A\, then 
K = [B, AB, A^B, . . .] has a full pivot structure J. 

(b) For every non-admissible full pivot structure J there exists an n x {m + n) matrix [B, A] having 
the full pivot structure J, for which K = [B, AB, A^B, . . .] does not have a full pivot structure. 

Proof, (a) Admissibility of J implies that B has a pivot- 1 column. Thus, the (infinite) controllability 
matrix K also has a pivot- 1 column, because i? is a sub-matrix of K. Now consider the induction 
hypothesis that the controllability matrix is known to have pivots at positions 1,2, ... ,k, with 1 < A; < 
n. From the admissible pivot structure of [B, A], either i? or ^ has a pivot-(A; + 1) column, depending 
on the value of jk+i- If Jfc+i < rn, then this column is in B hence it also appears in K. Otherwise, 
column jj^^i of [B, A] is in fact column Pk+i = jk+i — m of A. Equivalently, si = qm+e = k -\- 1 
for £ = pk+i- Because of the staircase structure of A (and because the prescribed pivot-1 column is 
in B) it holds that £ < k. Since K has a pivot-£ column according to the induction hypothesis, the 
matrix product AK now has a pivot-(A: + 1) column because of the staircase structure of A. But AK 
is a sub-matrix of K, whence it follows that K has a pivot-(A; + 1) column. This shows the induction 
step. Hence the controllability matrix has a full pivot structure. 

(b) See Appendix lAl □ 



(i) For an admissible pivot structure J for [B, A] there is a uniquely determined full pivot structure J 
which applies to every controllability matrix K that may occur for each arbitrary matrix [B, A] having 
the structure J. One can easily calculate J using the numbered Young diagram technique described in 
the following section. It is most clearly displayed in K for the example where each pivot-A; column in 
[B, A] is set equal to and each non-pivot column is set to zero. 

(ii) For given m and n, the total number of different admissible full pivot structures can be computed 
from the fact that an admissible pivot structure is completely determined by the pivot structure for B. 



Remarks. 
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4 The Young diagram associated with an admissible pivot structure 



Starting from an admissible pivot structure J for [B, A\ we now want to analyze the full pivot structure 
J for the (finite) controUabiUty matrix K = [B, AB, . . . , A^^^B] induced by J and describe their 
relation. 

Admissibility of J implies that 1 < ji < m, so that si is either zero (in which case A has no 
pivots) or si > 1. Together with the staircase form of A this means that for all = 1, . . . , n either 
Sk > k or Sk = 0. The sequence {k, S{k), S^{k), S^{k), . . .} therefore is strictly monotonically 
increasing until at some point the value occurs after which the sequence remains zero. This happens 
when S^{k) attains a value in {pA + 1, . . . , n}. Conversely, starting from a value i > the sequence 
{i, P{i), -P^(^), P^(^), . . .} is strictly monotonically decreasing until at some point the value occurs 
after which the sequence remains zero. This happens when P*(^) attains a value in {qi,... , q^n}. 

In this way, an admissible pivot structure J for [B,A] generates a uniquely specified full pivot 
structure J for the controllability matrix K. To visualize this, it is helpful to introduce an m x n 
array Y = {yij), defined as follows: entry y^j denotes the pivot position of vector i in the j-th block 
A^''^B of K (so that J{k) = {j — l)m + i where i and j are such that yij = k). In terms of the 
column-oriented description Q = {qi,q2,..., 5nm} of the pivot structure of K associated with the 
row-oriented full pivot structure J, it simply holds that yij = for alH = 1, . . . ,m and 

j = 1, . . . , n. The array Y can therefore be regarded as an m x n matrix representation of Q which 
allows a clearer expression of the role played by the block-partitioning of K. Obviously, there is a 
one-to-one correspondence between such an array Y (with entries in {0, 1, ... , n}) and the function Q 
(from {!,..., nm} to {0, 1, ... , n}). 

Example. Let m = 4, n = 6 and consider the admissible full pivot structure J = {3, 1, 5, 6, 4, 7} 
and its associated column-oriented description Q = {2,0,1, 5, 3, 4, 6, 0, 0, 0} for the 6 x 10 partitioned 
matrix [B, A] given by: 



[B,A] 



* 
+ 







+ 









* 

* 

+ 





* 
* 

* 

+ 






* 
* 

* 

+ 



Then the successor function S is given by S{0) = and {si, S2, S3, S4, S5, sq} = {3, 4, 6, 0, 0, 0} and 
the predecessor function P is given by -P(O) = and {pi,P2,P3,P4,P5,P6} = {0, 0, 1, 2, 0, 3}. Note 
that the matrix K =[B,AB, A^B, . . .] is of the form: 



K 



* 
+ 







+ 









* 

+ 





* 

* 

+ 






* 

* 

+ 







* * 



* 

* 
* 

+ 



This shows that the induced full pivot structure J for K is given by J = {3, 1, 7, 5, 4, 11} and it has 
an associated column-oriented description Q = {2, 0, 1, 5, 4, 0, 3, 0, 0, 0, 6, 0, . . . }. The corresponding 
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4x6 array Y is filled with the values in Q column after column, yielding the diagram: 



Y = 



2 


4 
































1 


3 


6 











5 


















Note that the first column of Y specifies the pivot structure of B, i.e.: {qi, q2, qs, q^} = {2, 0, 1, 5}. 
The other entries of Y satisfy the rule Uij+i = S{yij). 

Theorem 4.1 Let J be an admissible full pivot structure for the block-partitioned matrix \B, A\, with 
an associated column-oriented description Q = {qi, . . . , g^+n} <^nd the successor function S given 
by S = {si, . . . , Sn} = {qm+i-, ■ ■ ■ , Qm+n} and 5(0) = 0. Then J induces a full pivot structure J for 
the (finite) controllability matrix K = [B, AB, . . . , A^'^B] which is specified in terms of the m x n 
array Y associated with Q as follows. ■ 

(i) yi,i = Qifori = l,...,m; 

(ii) Vij+i = S{yij)for i = 1, . . . , m and j = 1, . . . , n - 1. 

Proof. As argued in the previous section, the admissible pivot structure J for [B, A] is entirely de- 
termined by the induced column-oriented pivot structure {qi, ... ,qm} for B. Given these (prescribed) 
pivot positions for the columns of B, the resulting pivot positions for the columns of the block AB 
are given by {S{qi), . . . , S{qm)}- Likewise, the pivot positions for the columns of the block A^B are 
given by {S'^{qi), . . . , S'^(gm)}. Proceeding in this fashion, it follows that the pivot structure J for 
K induced by J corresponds to an array Y which is described by: (i) the first column of Y, which 
corresponds to B and satisfies yj i = for i = 1, . . . , m; (ii) the other columns of Y, which are given 
by the recursion t/ij+i = S{yij) for i = 1, . . . , m and j = 1, . . . , n — 1. In part (a) of the proof of 
Theorem 13.41 it has already been argued that J obtained in this way describes a full pivot structure for 
K. □ 



The array Y in the theorem above has the property that the values 1, 2, . . . , n all occur precisely once 
while the other (m — l)n entries are all zero. The set of arrays Y with this property is denoted by 
y{m, n). Clearly, there is a one-to-one correspondence between this set of arrays and the set of full 
pivot structures for finite controllability matrices K of size n x nm. However, not all the arrays Y in 
the set y{m, n) are induced by some admissible pivot structure J for [B, A\. The following definition 
serves the goal of characterizing the subset of 3^(m, n) of arrays Y that are induced by admissible pivot 
structures. 

Definition 4.2 An array Y £ y{m,n) is called an admissible numbered Young diagram if it has the 
following three properties: 

( i) for all i = 1, . . . ,m and j = l,...,n — lit holds that yij+i > implies yij > ; 

(ii) the values n — pB -\- 1, ■ ■ ■ ,n all occur in different rows ofY as their last nonzero entries, where 
Pb is the number of nonzero rows ofY; 

(Hi) for all = 1,. . . ,m and = 1, . . . ,n — 1 it holds that yij+i > yi'j'+i > implies 

yi,j > Vi'j' > 0. 

Note that the number of nonzero rows of the array Y corresponding to the induced full pivot structure 
J in Theorem 14.11 is equal to the number of nonzero entries in the first column of Y, which is equal 
to the number of pivots in the matrix B. This explains the notation pB in the definition above. The 
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terminology 'numbered Young diagram' will become more clear below, when the relationship with 
nice selections and Young diagrams is explained. 

Theorem 4.3 (a) Let J be an admissible full pivot structure for the block-partitioned matrix [B, A\, 
then the induced full pivot structure J for the controllability matrix K = [B, AB, . . . , A^~^B] corre- 
sponds to an admissible numbered Young diagram Y. 

(b) Let Y be an admissible numbered Young diagram. Then there exists an admissible full pivot struc- 
ture J for [B, A] which induces the full pivot structure J for K = [B, AB, . . . , A^''^^B] which corre- 
sponds to Y. 

Proof, (a) From Theorem 14.11 we have that J induces the full pivot structure J for K which corre- 
sponds to an array Y G y{m,n) given by: (i) y^^i = for i = 1, . . . ,m; (ii) Uij+i = S{yij) for 
i = 1, . . . ,m and j = l,...,n — 1. Clearly, the z-th row of Y is entirely zero if and only if qi = 0. 
Hence the number of nonzero rows of Y is equal to the number of (prescribed) pivots in B. As we 
have seen, admissibility of J implies that Sfc = if and only if k ^ {pA -\- 1, . . . ,PA + Pb = n}. This 
shows that the last nonzero entries in the pB nonzero rows of Y have the values n — ps -\- 1, ■ ■ ■ ,n 
and they necessarily all occur in different rows. Next, if Vij+i > 0, then Vij+i = S{yij) with 
Uij > because S{0) = 0. This relationship is also described by the predecessor function P as 
Tjij = P(yj j+i) > 0. Note that in fact the restricted positive functions and P+ describe this 
relationship and they are both strictly monotonically increasing because of the staircase property of A. 
Therefore, by application of P+, the relationship Vij+i > Vi'j'+i > implies that yij > yj/ j/ > 0. 
This shows that Y is an admissible numbered Young diagram. 

(b) Suppose that Y € y{m, n) is an admissible numbered Young diagram. Consider the pB nonzero 
rows of Y. According to property (ii), the last nonzero entries of these rows precisely cover the range 
{n —pB + l, • • • , n}. It follows that all the other entries of 1" are < n —pB because every positive value 
from {1, . . . , n} occurs exactly once. Now consider the function S : {0,1, ... ,n} — > {0,1, ... ,n} 
defined from the values in Y as follows: S{0) = 0, S{yi,j) = yi.j+i for all i = 1, . . . ,m and 
j = I, . . . ,n — 1, and S{yi^n) = for alH = 1, . . . , m. Note that the pattern of positive values in the 
array Y is left-aligned according to property (i). This makes that the definition 5(0) = is consistent 
with the prescription S{yij) = yij+i in situations where yij = 0, and also with the prescription 
S{yi,n) = in situations where = 0. Note also that S{k) > for all = \, . . . ,n — pB and 
S{k) = for /c = 71 — p5 + 1, . . . , n (as well as for k = 0). The associated function is a bijection 
with domain {1, . . . , n — pb}- 

Property (iii) of Y now implies that 5+ is monotonically increasing. To see this, choose positive 
integers k and I with S{k) > S{i) > 0. Then choose the unique integers i, i', j and j' such that 
yi,j+i = S{k) and j/j'j'+i = S{£) and invoke property (iii) to obtain that k > £ > 0. Consequently, 

can be used to prescribe a staircase form for the matrix A. The positive values in {1, . . . , n} not 
occurring in the range of 5+ are precisely those occurring in the first column of Y. This first column 
of Y serves to describe a pivot structure for B. Together with this determines a full pivot structure 
J for [B, A] in which A has a staircase form. For J to be admissible, it remains to be shown that B has 
a prescribed pivot- 1 column, or equivalently that one of the entries in the first column of Y is equal to 
1. To see this, suppose that for some yij > it holds that S{yij) = yi,j+i = 1- Then the bijection 
5+ can only be monotonically increasing if S{1) = so that 1 does not belong to the domain of 5+, 
which requires 1 to belong to the set of pB largest values {n — + 1, . . . , n}. But then yi j > 1 also 
belongs to this set and occurs in a different row of Y, producing a contradiction. □ 

We thus have established a bijection between admissible pivot structures J for [B, A] and admissible 
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numbered Young diagrams Y associated with K. To relate these results to the well-known theory of 
nice selections and dynamical indices, the following definition is useful. 

Definition 4.4 Let m and n be given positive integers. 

(a) The set 'D{m, n) is defined as the set of all multi-indices d = [di ,d2,..., dm) G for which 
di + d2 + . . . + dm = n. 

(b) A selection of n columns from an n x nm controllability matrix K = [B,AB, . . . ,A^''^^B\ is 
called a nice selection if there exists a multi-index d G T>[m^ n) for which the selected set of columns 
is given by {A^~^Bei | j G {1, 2, . . . , di] for i = 1, 2, . . . , m}. 

(c) A nice pivot structure J for K is a full pivot structure for K which constitutes a nice selection of 
columns from K. 

(d) If J is a nice pivot structure for K, then the associated multi-index d G T){m, n) is called the vector 
of dynamical indices and each number di is called the i-th dynamical index (i = 1,2, ... ,m) of the 
nice pivot structure, or of the input pair {A, B). 

Nice selections and vectors of dynamical indices d are useful and well-known concepts for studying 
the rank structures that can be exhibited by a controllability matrix K. The most well-known nice 
selection is the Kronecker nice selection, which consists of the first n linearly independent columns of 
K. Every nice selection may occur as the Kronecker nice selection for some controllability matrix K. 
(Cf., e.g., Hi and the references given there.) In the concept of nice selections though, there are no a 
priori rank or linear independence requirements and no triangularity conditions. Conversely, for a nice 
pivot structure it is not required that the column selection is a Kronecker nice selection. Note also that 
there are n! different nice pivot structures all corresponding to the same nice selection. 

Above it has been shown that an admissible pivot structure for [-B, A] induces a corresponding 
full pivot structure for K for which the associated array Y G 3^(m, n) is an admissible numbered 
Young diagram. Conversely, all admissible numbered Young diagrams are induced in this way. An 
admissible numbered Young diagram specifies a selection of n columns of K, which constitutes an 
upper triangular sub-matrix; therefore these n columns are linearly independent. From the definition 
of a nice selection it should be clear that any nice selection can be represented by an m x n binary array 
Z = (zij) in the following way: Zij = 1 if column i of the j-th block A^^^B of K is included in the 
nice selection, and Zij = otherwise. The nonzero entries in such an array Z exhibit a left-aligned 
pattern and the dynamical index di denotes the number of nonzero entries in the i-th row of Z, while 
di -\- . . . -\- dm = n. Such an array Z is closely related to the concept of a Young diagram, see [2 |. As 
we have seen, any admissible numbered Young diagram Y is left-aligned and it therefore gives rise to 
an associated nice selection; the induced full pivot structure J is a nice pivot structure for K. This also 
explains our terminology. For the purpose of the design of local canonical forms for various classes of 
linear multivariable systems, it is important that there exists an admissible numbered Young diagram 
for every nice selection. We therefore continue to study the relationship between nice selections and 
admissible numbered Young diagrams. 

Let Z be a Young diagram, i.e., a left-aligned m x n binary array corresponding to a nice selec- 
tion with an associated vector of dynamical indices d = {di, . . . , dm)- A numbered Young diagram 
is obtained from Z by replacing the unit entries in Z by the numbers 1, 2, . . . , n in some arbitrary 
order, so that they all occur exactly once. The set of m x n numbered Young diagrams is the subset of 
y{m, n) of left-aligned arrays. We will now show that for every Young diagram Z there exists an asso- 
ciated admissible numbered Young diagram Y. More precisely, we will characterize all the admissible 
numbered Young diagrams Y that correspond to Z. 
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To do this, it is convenient to associate with every left-aligned array Y G y{m, n) a corresponding 
right-aligned array Yr as follows. If Y is left-aligned then this means that there is an associated vector 
of dynamical indices d = {di, . . . ,dm) such that yij > iff j < dj. Thus, the i-th row of Y has 
positive entries at its first di positions and zero entries at the remaining n — di positions. Then Yj. is 
defined by: {Yr)ij = for 1 < j < n — and {Yj.)ij = yij-n+d, for n — dj + 1 < j < n. In other 
words: the di positive entries in the i-th row are all shifted n — di positions to the right. 

Proposition 4.5 Let Z be an m x n Young diagram corresponding to a nice selection with an asso- 
ciated vector of dynamical indices d = (di, . . . , dm)- An m x n left-aligned array Y corresponding 
to the same vector of dynamical indices d, is an admissible numbered Young diagram if and only if 
there exists an m x m permutation matrix 11 for the associated right-aligned array Yr such that the 
nm-vector vecilVYr) = ((ny^ei)^, (111^^62)^, • • • , (ni^e„)-^) G R""* obtained by stacking the 
n columns of the array TYYr, has the property that if the zero entries are deleted then the n-vector 
(1, 2, 3, ... , n)^ is obtained. 

Proof. Suppose Y is an admissible numbered Young diagram corresponding to the vector of dynam- 
ical indices d. Consider the ps = n — pA nonzero values in the last column of Yr (where pB denotes 
the number of pivots in B, i.e. the number of nonzero entries in {qi, ... ,qjn} which also is the number 
of nonzero rows in Y as well as in Yr). These values constitute a permutation of the set of values 
{pa + 1, . . . ,pa + Pb = n}. Now consider the predecessors {P{pA + l)i • • • , Note that the 

nonzero values among these predecessors show up in an increasing order, because P+ is monotoni- 
cally increasing. Repeating the argument, it follows that the same permutation of the nonzero rows 
of Yr which makes that the nonzero entries in its last column appear in an increasing order, achieves 
that such a property holds for each of the columns of Yr. Consequently, when all the columns of the 
row -permuted array Yr are stacked into a vector with nm entries using the well-known vec(-) operator, 
a column vector remains which is equal to (1, 2, ... , n)'^ when all the zeros entries are deleted. 

Conversely, starting from the given vector of dynamical indices d and an arbitrary choice of 11 
permuting the nonzero rows of Z, the nm-vector with the given property and the right-aligned arrays 
HYr and Yr and the left-aligned array Y are completely determined. The left-alignment property (i) of 
an admissible numbered Young diagram Y is built-in. Properties (ii) and (iii) of Y are not difficult to 
verify either, because they are easy for 111^ and Yr and shifting the rows to move between Yr and Y 
does not basically change the requirements (one only needs to take into account that zeros may occur 
to the left of a string of nonzero entries in Yr, but the dynamical indices now specify the length of such 
a string in advance). □ 

Note that the technique used in the proof of this proposition is constructive and can be used to generate 
all the admissible numbered Young diagrams corresponding to a given nice selection. There are ps! 
different possibilities, where pb can be read off from d as the number of dynamical indices di > 0. 

For given m and n and for each nice selection with a vector of dynamical indices d, one can consider 
the family J^{d) of controllable pairs which have the additional property that the selected columns from 
the controllability matrix are linearly independent. Then we know that each controllable pair {A, B) 
Ues in at least one of the families F{d), d e'D{m,n). 

Now consider the family of all controllable pairs (A, B) with A discrete-time asymptotically stable, 
and the question of how to parameterize this family up to state isomorphism. (I.e., up to multiplication 
of the controllability matrix if by a nonsingular n x n matrix on the left.) Every such pair [A, B) 
corresponds to a positive definite controllability Gramian Wc, which can be factored into Wc = M'^M 
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by making a well-defined choice for M, e.g. by prescribing it to be a positive upper triangular Cholesky 
factor. Such a choice can be parameterized to involve precisely n(n+l)/2 independent real parameters. 
Using M to act as a state isomorphism transforms {A, B) into an input-normal controllable pair, but it 
does not affect any linear dependence relations between the columns of the controllability matrix K. 
Hence it allows one to restrict attention to the question of parameterizing the family of input-normal 
controllable pairs {A, B) up to orthogonal state isomorphism. 

Note that an input-normal controllable pair {A,B) corresponds to a row-orthonormal partitioned 
matrix [-B,^] for which A is asymptotically stable, and vice versa (see e.g. Q, ||8l). Then for each 
admissible numbered Young diagram Y the family of row-orthonormal [B, A\ with A asymptotically 
stable and with an admissible pivot structure corresponding to Y , forms a local canonical form for 
this family. This set of local canonical forms is covering this family in the sense that for each row- 
orthonormal [B, A] with A asymptotically stable there exists an admissible numbered Young diagram 
Y and an orthogonal matrix Q such that [QB, QAQ^] has the admissible pivot structure associated 
with Y. Furthermore, because of uniqueness of the associated QR-decomposition due to positive upper 
triangularity, for such a combination of [B, A\ and Y the matrix Q and hence [QB, QAQ^] is unique. 

An interesting question is how to obtain a minimal sub-atlas of this atlas of local canonical forms, 
minimal in the sense that no further local canonical forms can be left out without losing the property 
of covering the family. To obtain a minimal sub-atlas we have to choose one of the local canonical 
forms for each d € P(m, n). This implies that for each d G D(m, n) we have to choose one of the 
Pb \ possible numberings of the associated Young diagram. As each such numbering is associated with 
a permutation of the nonzero rows of the Young diagram this choice can be fixed by specifying that 
permutation. One possible choice is the unique permutation for which the permuted dynamical indices 
form a non-increasing sequence, while the order of the rows which have the same dynamical index 
is kept the same. Note that this permutation is used only to determine the numbering in the Young 
diagram, the ordering of the dynamical indices is left unchanged. With hindsight one can say that this 
particular choice to obtain a minimal atlas was used in [61 in a similar approach for continuous-time 
input-normal pairs and lossless systems. Just as in that paper for the continuous-time case, here each 
local canonical form on discrete-time asymptotically stable input normal systems defines a balanced 
local canonical form on minimal discrete-time lossless systems of order n. How these balanced local 
canonical forms for minimal discrete-time lossless systems of order n are related to those constructed 
in 171 by means of the tangential Schur algorithm is the topic of the next section. 

5 Atlases of balanced canonical forms for lossless systems 

We now have two approaches to arrive at an atlas of overlapping balanced canonical forms for discrete- 
time lossless systems: one using the balanced realizations associated with the tangential Schur algo- 
rithm and one based on balanced realizations with an imposed pivot structure on the row-orthonormal 
matrix [i?, A], hence on the orthogonal realization matrix R. However, one of our main results is that 
the second approach corresponds to making special choices for the direction vectors in the first ap- 
proach. Hence the atlas of overlapping balanced canonical forms resulting from the second approach 
is a sub-atlas of the atlas of overlapping balanced canonical forms in the first approach. The precise 
formulation is as follows. 

Theorem 5.1 Let Y be an admissible numbered Young diagram, corresponding to an associated nice 
pivot structure J (for controllability matrices) and an admissible pivot structure J (for n x [m + n) 
matrices). For each k = 1,2, ... ,n, choose the direction vector Un+i~k equal to ej(^.), the i{k)-th 
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standard basis vector in M."^, where {i{k), j{k)) denotes the unique pair of indices such that y^i^'jj(^f.-j = 
k. Then for any choice of the Schur parameter vectors vi,V2, ■ ■ ■ ,Vn (all of length < 1) and for any 
choice of the orthogonal matrix Dq, consider the [m + n) x (m + n) orthogonal realization matrix R 
given by (|£l). It follows that J is an admissible pivot structure for the sub-matrix [B, A\ and J is a nice 
pivot structure for the controllability matrix K. 

A detailed technical proof of this theorem is given in Appendix |B] 

From the point of view of the tangential Schur algorithm, it is of interest also to directly characterize 
all the sequences of direction vectors ui,U2, ■ ■ ■ ,Un that give rise to an admissible pivot structure for 
the matrix [B, A] (and an accompanying nice pivot structure for the controllability matrix K). 

Theorem 5.2 Consider a chart associated with the tangential Schur algorithm ( with all the interpola- 
tion points Wk located at the origin), specified by a sequence of direction vectors {ui, . . . , Then 
each [B,A\ resulting from this chart exhibits an admissible pivot structure, irrespective of the choice 
of Schur vectors vi, . . . ,Vn, if the sequence of direction vectors consists of standard basis vectors, say 
Uk = &fi{k) for some indices . . . , ^^{n) chosen from {1,2,..., m}, satisfying the following con- 
dition: 

for each k = 1,2,... ,n — 1, if there exists a largest index I strictly less than k such that Ui = e^(,fc), 
then ^{k + 1) is from the set {^{i + 1), . . . , /^(A;)}. 

Proof. This follows directly from the properties of the three procedures introduced in Appendix iBl to 
generate an admissible numbered Young diagram. Details are left to the reader. □ 

Example. Consider the same situation as for the example in Appendix |Bj where m = 5, n = 12 and 

(^i) 92, 93j 94) 95) = (4, 1, 9, 0, 7). There it is remarked that the choice of direction vectors Un+i-k = 
ei(^k) can be rewritten as Uk = e^{k) where n{k) denotes the index of the value 1 in the vector rjk, 
generated by the 'third procedure'. In this example it follows that the sequence {ui,U2, ■ ■ ■ , U12} is 
given by {e2, 65, ei, 63, 62, 65, ei, 62, ei, 62, 62, 62}- Note that this sequence satisfies the condition of 
Theorem 15.21 for all A; = 1,2,..., 11. E.g., for = 6 the previous occurrence of the vector uq = 65 
happened for i = 2. The condition of the theorem requires U7 to occur in the set {us, U4, U5, ug} = 
{ei, 63, 62, 65}, which indeed is the case. 

6 Conclusions and discussion 

In this paper we have developed a detailed procedure to construct an atlas of overlapping (local) bal- 
anced canonical forms for MIMO discrete-time lossless systems {A, B,C,D) of a given order n. To 
achieve this, the concept of an admissible pivot structure for [B, A] has been introduced, which induces 
a full (nice) pivot structure for the associated controllability matrix K. The approach taken in Sections 
[3]and|4]has a wider range of applicability though: it builds on the algebraic relationship between [B, A] 
and K, and it neither relies on input-normality of {A, B) nor on the discrete-time setting for the loss- 
less systems. When one is dealing with a system having a special (non-generic) Kronecker structure, 
this can be recognized in these local canonical forms by certain entries becoming zero. To demonstrate 
the structure of the charts that constitute the atlases discussed in this paper, a detailed worked example 
is included in Appendix ICl 

One of the main practical reasons for studying overlapping canonical forms is that (iterative) iden- 
tification and optimization algorithms in the MIMO case may suffer from numerical ill-conditioning 
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and slow convergence when they pass by systems that not too far from systems with a non-generic 
structure. Switching charts may then help to improve the algorithmic performance. The connection 
in Section [5] with the atlas of charts developed for discrete-time input-normal pairs involving the tan- 
gential Schur algorithm is useful, because that set-up involves (well conditioned!) orthogonal matrix 
computations while it is tailored to deal with the important class of stable systems. The tangential 
Schur algorithm provides one with a lot of flexibility to design local canonical forms. Since it is com- 
putationally expensive to switch charts at each and every iteration, a suitably chosen finite sub-atlas 
is welcome. In the present paper we have indicated the restrictions that should be taken into account 
when choosing direction vectors Uk from the set of standard basis vectors, if a pivot structure is to show 
up not only in [B, A] but also in the controllability matrix K. When a nice pivot structure is present in 
K, this has the advantage that the impact of state vector truncation is easier to analyze; controllability 
is then preserved. This is of importance in the context of model order reduction applications. 

Future research addresses the issue of monitoring the conditioning of a chart (i.e., a local canonical 
form) at a given system, and the issue of selecting a better chart when switching becomes necessary. 
Since the total number of charts in an atlas quickly grows large with the dimensions m and n (even for 
the case of admissible pivot structures) it may not be attractive to carry out a full search for a better 
chart over the entire atlas. The rank structure in K can then be instrumental in designing a quick on- 
line algorithm which guarantees a certain degree of conditioning improvement. This is currently under 
investigation. 
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A Proof of part (b) of Theorem 1X4 

Consider a non-admissible full pivot structure J = {j'l, j2, ■ ■ ■ ,jn} for [-B,^]. Then it either holds 
that: (1) J does not prescribe B to have have a pivot- 1 column, or (2) J does prescribe B to have a 
pivot- 1 column, but J does not impose a staircase structure on A. 

In case (1) it holds that ji > m. We distinguish between two situations, (i) If ji = m + 1, then the 
first column of yl is a pivot-1 column. Then consider the following example: for each /c = 1, . . . , n 
let column of [B, A] be equal to the fc-th standard basis vector G R" and let the remaining m 
columns of [B, A] all be zero. Clearly, [B, A] exhibits the given full pivot structure J, but ejB = 
and ejA = ef so that ejK = 0. Hence K does not have a pivot-1 column, so it does not have a full 
pivot structure. 

(ii) If instead ji > m + 1, then consider the following example: for each k = \, . . . ,n choose column 
jk of [B, A\ to be a positive pivot- A; vector and let the remaining m columns of [B, A\ all be chosen 
to be strictly positive (so that effectively they are all positive pivot-n vectors). Clearly, [B, A] exhibits 
the given full pivot structure J. Note that each column in B is (effectively) a positive pivot-fc vector 
with k > 2. Now, if u is a positive pivot-/c vector, then Av is a positively weighted linear combination 
of the first k columns of A. Since all columns of A are (effectively) positive pivot-i' vectors for certain 
values of i, the vector Av is a positive pivot-p vector where p is the maximal (effective) pivot position 
among the first k columns of A. Now, the first column of A has at least two nonzero entries, because 
ji > m + 1. Therefore, each column of AB is (effectively) a positive pivot-p vector with p > 2. 
By induction it follows that all columns of K are (effectively) positive pivot-p vectors with p > 2. 
Consequently, K does not have a pivot-1 column, so it does not have a full pivot structure. 
In case (2) it holds that ji < m, but the staircase structure does not necessarily hold for A. We again 
distinguish between two situations, (i) Suppose that for some k < n there is a pivot-A; vector in A for 
which there is either a non-pivot column in A preceding it, or a pivot-^ vector preceding it with £ > k. 
Then consider basically the same example as used in case (1) part (ii): for each k = 1, . . . , n choose 
column jfc of [B, A] to be a positive pivot-/c vector and let the remaining m columns of [B, A] all be 
chosen to be strictly positive (so that effectively they are all positive pivot-n vectors). For this example 
it now follows that K does not have a pivot-A; vector, because B does not have one and because for 
all p = 1, . . . , n the maximum (effective) pivot position among the first p columns of A can never be 
equal to k. 

(ii) For all A; < n the pivot-A; vectors in A respect the staircase structure, but there is a prescribed 
pivot-n vector in A which is directly preceded by a non-pivot column. If this pivot-n vector occurs in 
the last column of A, then one may consider the same kind of example as used in case (1) part (i): for 
each k = 1, . . . , n let column of [B, A] be equal to and let the remaining m columns of [B, A] 
all be zero. Now eJ^B = and e^A = so that eJ^K = 0. Hence K does not have a pivot-n column, 
so it does not have a full pivot structure. If the pivot-n vector does not occur in the last column of A, 
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then the last column of ^ is a non-pivot column. Summarizing, we then are in a situation where the 
first PA — i columns of A exhibit a staircase structure, column p := jn — m > pA of A is a. pivot-n 
column and the two columns p—1 and nof A are both non-pivot columns. Then consider the following 
example. Column jp_i of [B, A] is defined as the pivot- (p — 1) vector |ep_i. Column jp of [B, A] is 
defined as the pivot-p vector fcp + ^Cp-i. The 3x3 sub-matrix S of A constituted by the intersection 

/ ^ 125 125 \ 

of its rows and columns with indices p — 1, p and n is defined as 5" = — i| it . All 

I _3 16 12 / 

\ 5 25 25 / 

remaining entries at pivot positions in [B, A] are defined to be 1 and all other entries are set to zero. 
For this example it will be shown that the entries in the last row of K are never positive, so that K 
does not have a full pivot structure. Note that all columns in B have a last entry that is equal to zero. 
All other columns in K are of the form Av for some vector G M". For Av to have a nonzero last 
entry, at least one of the entries in positions p — l,p and n of u must be nonzero. Such vectors v must 
come from repeated pre-multiplication of the columns of B by the matrix A. The first vectors to have 
such a structure are the pivot- — 1) vector and the pivot-p vector that both occur among the columns 
of B and the first pA — 1 columns of A. Once such vectors v are multiplied by A, only the entries in 
positions p — l,p and n can become nonzero: the subspace spanned by Sp-i, Cp and e„ is an invariant 
subspace of A. Restricting to this subspace, the matrix A is represented by the sub-matrix S given 
above. Consequently, the entries in the last row of K are either zero or obtained as the entries in the 

/ 3 16 \ 
/ 5 25 \ 

last row of the controllability matrix of the pair (T, S) with T = | which represents the 

V / 

pivot- (p — 1) and the pivot-p vector in this new notation. For the matrix S it is easily established that 
S = It therefore suffices to compute the last row of the matrix [T, ST, S'^T], which is equal to 

(0, 0, — 0, — ~^)- This proves that all these entries are indeed non-positive. □ 

Remark. 

In the case of balanced realizations of lossless systems we will in addition require [S, A] to have 
orthonormal rows. The proof above does not entirely apply to this restricted situation. For example if 
[B, A] only has non-negative entries then orthogonality of the rows requires that in each column there 
is at most one nonzero entry. This requirement is violated by the counterexamples presented in case 
(1) part (ii) and in case (2) part (i) of the proof, because non-pivot columns are chosen to be positive 
pivot-n vectors. How to obtain a proof for this more restricted orthonormal case is an open problem at 
this point. Note however that the counterexamples presented in case (1) part (i) and in case (2) part (ii) 
of the proof have in fact already been designed to involve [B, A] with orthonormal rows. 



B Proof of Theorem O 

To prove this theorem it is helpful first to reconsider the precise relationship between Y, J, J and the 
pivot structure {qi, ... ,qm} of the matrix B. Recall that the number 1 appears in the pivot structure 
{qi, ... ,qm} for B because of admissibility and this sequence completely characterizes the successor 
function S, the column-oriented description Q, the admissible pivot structure J, the nice pivot structure 
J and the admissible numbered Young diagram Y in the way explained before. 



We have previously introduced the following construction procedure for Y from {qi, . . . , q-m}'- 
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Procedure 1 

(a) Construction: the successor function S is defined as the increasing sequence of all positive num- 
bers in {1,2, ... ,n} not occurring in {gi, . . . , g^} completed by a sequence of pB zeros; in addition 
S{0) = 0. 

(b) Initialization: set y^^i := for i = 1, . . . , m. 

(c) Recursion: set yij+i := S{yij) for i = 1, . . . , m and j = 1, . . . , n — 1. 

A second way to generate Y in a dynamical fashion which avoids the exphcit construction of S, is by 
means of the following procedure: 
Procedure 2 

(a) Initialization: set j/j^i = qi for i = 1, . . . , m. 

(b) Recursion: for A; = 1, 2, . . . , n, if the value k has not yet been assigned to an entry of Y then select 
the smallest nonzero number in Y for which the entry j/jj+i immediately to its right has still not 

been assigned some value and set Vij+i = k. 

(c) Termination: set all the remaining entries of Y equal to zero. 

It is not very hard to establish that the array Y constructed in this fashion is indeed admissible and 
identical to the array Y constructed previously with the help of S. 

A third way to generate Y from {gi, . . . , q^} employs a sequence of vectors r]k {k = f),l,2, . . . , n) 
and proceeds as follows in a backward fashion: 
Procedure 3 

(a) Initialization: set r)n = {qi, ■ ■ ■ , Qm)'^ ■ 

(b) Backward recursion: for k = n — 1, ... ,1,0 construct r]k from ij^+i by executing the following 
three rules in the given order: 

(1) if {Vk+i)i = then set (r/fc), := 0; 

(2) if {rik+i)i > 1 then set {rik)i := {vk+i)i - 1; 

(3) if {i]k+i)i = 1 then define ^j. as the smallest positive number different from all the entries of r]k 
already assigned by rules (1) and (2); if ^ k then set {r)k)i := else set {rik)i := 0. 

(c) Construction: for each i = 1,... ,m consider the di values of k for which {'r]n+i-k)i = 1 and 
assign these values (in increasing order) to the first di entries of row i of Y; set all other entries to 0. 

The validity of this third procedure for generating Y can be seen as follows. First, note that because of 
rule (2) in each recursion step (b), the first column of Y attains the required form containing qi,. . . , qm< 
since the number 1 first occurs in position i of rjn+i-k for k = qi. Also note that the positive integers 
in r]n are all different (since this holds for gi, . . . , q„i) and that the rules in each recursion step (b) are 
such that this property is preserved for all vectors rik- Next, these rules are such that each vector rjk 
has precisely one entry equal to 1, for all A; = 1, . . . , n. The construction in step (c) is such that all the 
numbers from {1, . . . , n} show up precisely once in a corresponding left-aligned numbered Young dia- 
gram Y. Finally, rule (3) in each recursion step (b) guarantees that yij is followed by yi,j+i = S{yij): 
note that yij = kis equivalent to {r]n+i-k)i = 1 and {r]n-k)i = Ck impUes that j/jj+i = k + ^k', here 
^jfc > is chosen as small as possible, precisely in Une with the second procedure for generating Y. 

Example. Consider the situation with m = 5, n = 12 and (gi, q2, qs, q^, gs) = (4, 1, 9, 0, 7). Then 
the successor function S is described by {si,S2, ■■■ , sio} = {2, 3, 5, 6, 8, 10, 11, 12, 0, 0, 0, 0}. It 
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follows that = 4 and the corresponding admissible numbered Young diagram Y is given by: 
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11 























where all the zeros are omitted for clarity. The second procedure for generating Y without the explicit 
construction of the successor function S yields the same result. It proceeds from the given first column 
of Y by putting the value 2 after the value 1, then the value 3 after the value 2, then the value 5 after 
the value 3, then the value 6 after the value 4, and so on. The values 1, 4, 7 and 9 are skipped, because 
they have already been assigned to the first column of Y. 

The third procedure for generating Y involves the backward recursion for the construction of the 
vectors r/^, for A; = n, n — 1, . . . , 1, 0. This produces the following sequence: 
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For instance, the vector r]i is obtained from the vector i]s as follows. First all the entries equal to zero 
are copied and all the values (r/g)^ > 2 are decreased by 1 to produce the corresponding values of {7]j)i. 
The value of (777)2 is addressed last, because (778)2 = 1- At that stage the values 1 and 2 have already 
been assigned to some entries of 777 and it holds that ^7 = 3. Because ^7 = 3 is not larger than the 
index k = 7, this value is assigned to (777)2. 

Once the vectors 7/12, t/h, . . . , 771 have been constructed, the array Y is constructed by considering 
the positions of the entries 1. For the first row, these positions are subsequently 4, 6 and 10 (proceeding 
in the given order from 7/12 to 771 corresponding to the index n + 1 — k). For the second row we have: 
1, 2, 3, 5, 8 and 12, and so on. 

The third way of characterizing Y in terms of the pivot structure {qi, . . . , qm} has a number of prop- 
erties that are worth noting in view of the proof of Theorem 15 . 1 1 below. First, note that rule (3) in step 
(b) implies that the maximum value among the entries of 77^ is at most k. (Therefore, 770 is the zero 
vector.) Second, all the positive entries of a vector r]^ are different. This makes that if (?7fe+i)i = 1 
then a positive value is assigned to (77^)4 for A; > pb and the value is assigned for A: < pb- Third, 
note that the sequence of values ^n-i,Cn-2, ■ ■ ■ t£,pb (^^ that backward order) is increasing. Fourth, 
the choice of direction vectors n„+i_fc = ej(fc) can be rewritten as Uk = efi{k) where /i(/c) denotes 
the index of the value 1 in the vector 7/^. Note that according to this notation, (7/^)^(^+1) = for 
k = n — l,n — 2, . . . ,pb- Finally, it will be shown that the vectors 77^ represent the pivot structures 
for the sequence of lossless systems of orders k = 1,2, ... ,n encountered in the tangential Schur 
algorithm for the particular choice of direction vectors specified in Theorem 15. II 

Proof of Theorem 15. II Consider the matrix product 

R = Tn---r^RoAj---Al. 
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Note that the product r„ • • • rii?o is positive m-upper Hessenberg for any choice of Schur vectors 
vi,. . . ,Vn- Post-multiphcation by the matrix only affects the last m + 1 columns, because the 
matrix is given by 





Af = 



In-l 











°M(1) 



where jLt(l) denotes the location of the entry 1 in the vector r/i which features in the third method for 
the construction of Y from {gi, . . . , g^}. The precise effect is as follows: 

(i) column n of r„ • • • Fii^o (having a pivot in its last position) is moved into column n + /x(l) — 1; 

(ii) column n + of • • • Fii^o is moved into column n + m; 

(iii) columns n + 1, . . . , n + jLt(l) — 1 and n + + 1, . . . , tt, + m of r„ ■ ■ ■ ViRq are moved one 
position to the left, into columns n, . . . , n + — 2 and n + . . . . n + m — 1, respectively. 
Note that the last row of r„ • • • Fii^oAj^ can be regarded to have the structure: 

[0 Si ^1 ] 

with Ai of size 1x1 and Bi of size 1 x m. The 1 x (m+1) partitioned matrix [Bi , Ai] has an admissible 
pivot structure for which the column-oriented pivot structure of Bi is given by {0, . . . , 0, 1, 0, . . . , 0} 
with the value 1 in position In other words, the pivot structure of Bi is described by 7]i. 

Consider the last k rows of the matrix product r„ • • • Fii^oA^ • • • A^. Note that these can be 
regarded to constitute the structure: 

[0 Bk Ak] 

with Af; of size k x k and Bj^ of size k x m. Now suppose that the k x {m + k) partitioned matrix 
[Bk, Ak] is known to have an admissible pivot structure for which the column-oriented pivot structure 
of Bk is given by the vector rjk- (This is the induction hypothesis.) We consider what happens under 
post-multiplication by the matrix A^^^. Note that this matrix is given by: 



In-k-l 










T 
















h 



Therefore, post-multiplication by A^_j^^ only acts on the columns n — A;, . . . , n — A; -|- m of the matrix 

F„---Fii?flAf •••A^. 

The partitioned matrix [B^+i, A^+i] is then formed as 



[Bk+ii^k+i] 



^ (3 a 
Bk Ak 











Im - e^(fc+i)e^(jt+i) e^(fe+i) 
4 



where 7 is a positive scalar and /3 and a are 1 x m and 1 x A; row vectors, respectively. It follows that 
the post-multiplying matrix carries out the following action: 
(i) the columns involving Ak remain unchanged; 
7 /3 a 



(ii) column 1 of 



Bk Ak 



having a pivot in its first position, is moved into column i^(k + 1); 
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7 




/3 



a 



(iii) column /x(A; + 1) + 1 of 

(iv) columns 2, . . . , ^{k + 1) and ij,{k + 1) + 2, . . . , m + 1 of 



is moved into column m + 1; 



7 




are moved one position 



/3 a 
Bk A.k 

to the left, into columns 1, . . . , ^{k + 1) — 1 and ^{k + 1) + 1, . . . , m, respectively. 
This shows that the pivot structure of is obtained from the pivot structure of Bk in the following 
way: all the nonzero entries of the structure vector -qk are increased by 1 except for the entry with index 
li{k + 1), which is reset to 1. This means that the pivot structure of Bk+i is indeed given by the vector 

rik+i- 

It remains to show that the matrix Ak+i again has a staircase form. Now, Ak has a staircase form 
according to the induction hypothesis and Ak+i is recognized to be of the form 



A 



k+l 



5 a 
e Ak 



where 5 is the /u(A; + l)-st entry of the row vector /3 and e is the ^{k + l)-st column of Bk- This 
means that the pivot in the first column of Ak+i shows up in position (rjk) fj,(k+i) + 1> which is equal to 
^k + 1- However, it has already been established that the sequence ^n-i,Cn-2, ■ ■ ■ , Cpb increasing. 
Therefore, Ak+i also has a staircase form. 

By induction this shows for all /c = 1, 2, . . . , n, that [Bk, Ak] has an admissible pivot structure for 
which the vector rjk specifies the pivot structure of the matrix Bk- In particular, for /c = n the claim of 
the theorem follows. □ 



C An atlas for input-normal pairs {A, B) under orthogonal state-space 
equivalence, with m = 3 and n = 4 

To illustrate the results and constructions of this paper, we here present an atlas for the manifold of 
(controllable) input-normal pairs {A, B) under orthogonal state-space equivalence, for the non-trivial 
case m = 3 and n = 4. Each of the charts in this atlas gives rise to a particular full pivot structure in 
the controllability matrix K and an admissible pivot structure for the row-orthonormal matrix [-B, A]. 
For given m and n, the number of different admissible numbered Young diagrams (see the end 

of Section (3]) is specified by ^^"^i™-"} ^'^'^^^^^ l) '^^ ^^^^ m = 3 and n = 4 this 
amounts to 39. To obtain a minimal sub-atlas, precisely one chart should be included for each nice 
selection, i.e. for each vector of dynamical indices d in 'D{m, n). The cardinality of D(?7i, n) is easily 

(Tfi -\- Tt — 1 \ 
^ j . For the case m = 3 and n = 4 this implies that a minimal sub-atlas 

consists of 15 charts. In Tables [T][3] the 15 different vectors of dynamical indices for this example are 
displayed, along with the corresponding 39 admissible numbered Young diagrams and their associated 
pivot structures in K and in [-B, A]. 

To arrive at an explicit parameterization of a chart in these tables, one may proceed in the discrete- 
time case by exploiting Eqn. ^ for the construction of orthogonal realization matrices, correspond- 
ing to balanced realizations of discrete-time lossless systems. Here the sequence of direction vectors 
{ui, U2, U3, U4} is chosen to consist of particular standard basis vectors, as indicated for each chart in 
these tables too. The parameters are then provided by the sequence of Schur vectors {wi, f2, fa, f4} 
which are all required to be of length < 1. The 3x3 orthogonal matrix block Dq can be set to 



22 



any fixed value; the choice Dq = /s is a convenient one. The latter is a consequence of the general 

D C 



fact that if an orthogonal realization matrix R 



B A 



is generated by Eqn. ^ for some Dq, 



yield the realization matrix 



, which exhibits exactly the same input pair (^4, B). 



{ui, . . . , Un} and {vi, . . . , Vn}, then the alternative choices Im, {ui, . . . , u„} and {DqVi, . . . , DqVu} 

DoD DoC 
B A 

From such a (minimal or non-minimal) atlas for input-normal pairs under orthogonal state-space 
equivalence, a corresponding atlas for all input-normal pairs of the given dimensions m and n is directly 
obtained by regarding the associated manifold as a Cartesian product of the previous manifold with the 
orthogonal group 0{n), related to the choice of state-space transformation. 

To arrive at a corresponding atlas for m x m lossless systems of order n one may instead regard this 
space as a Cartesian product of the previous manifold with the orthogonal group 0{m), now related to 
the choice of Dq. 

To obtain an atlas for asymptotically stable discrete-time systems of order n with m inputs and p 
outputs, one may proceed by taking all the entries of C and D (of sizes p x n and p xm, respectively) 
to be free parameters, only subject to the constraint that observability needs to hold for the pair (C, A) 
(a property which is then generically satisfied in each chart, i.e. it only excludes a thin subset of 
parameter vectors). Such an approach is useful in system identification, for instance in conjunction 
with the method of separable least-squares (see [1]). Then we may have to consider output-normal 
forms instead, but this can be achieved easily using input-output duality. 

Finally, to deal with the continuous-time case, the well-known bilinear transform can of course be 
applied. However, this will in general destroy the pivot structure in K and in [-B, To employ the 
results directly in the continuous-time case too, note that the pivot structures for (controllable) input- 
normal pairs [B, A] as given in the Tables [Tl[3] do in fact apply to the continuous-time case already, 
giving rise to local canonical forms that can be computed numerically for a given state-space realization 
in a straightforward way. What at present seems to be lacking in the continuous-time case is an explicit 
parameterization of these local canonical forms (such as may be required in system identification). 
This is currently the topic of ongoing research. 
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Chart Young diagram Z and Admissible nmnbered Full pivot Sequence of direction Admissible pivot 
dynamical index d Young diagrams Y structure 3 for K vectors {t^i , U2, 'us, ^^4} structure J for [B, A] 



Structure of 
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{1,4,7,10} {ei,ei,ei,ei} {1,4,5, 
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(4, 0, 0) 



{1,4,7,2} {e2,ei,ei,ei} {1,4,6,2} 









+ * 



(3,1,0) 



{1,4,2,7} {ei,e2,ei,ei} {1,4,2,5} 



+ * 



* * * 
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{1,4,7,3} {e3,ei,ei,ei} {1,4,6,3} 
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(3,0,1) 



{1,4,3,7} {ei,e3,ei,ei} {1,4,3,5} 
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{1,2,4,5} {e2,ei,e2,ei} {1,2,4,5} 



(2,2,0) 
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{2,1,5,4} {61,62,61,62} {2,1,4,5} 
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{1,4,2,3} {63,62,61,61} {1,4,2,3} 



{1,4,3,2} {62,63,61,61} {1,4,3,2} 



{1,2,4,3} {63,61,62,61} {1,2,4,3} 



{1,3,4,2} {62,61,63,61} {1,3,4,2} 



{1,2,3,4} {61,63,62,61} {1,2,3,4} 



{1,3,2,4} {61,62,63,61} {1,3,2,4} 
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Table 1: Charts 1-5 for [-B, for the case m = 3 and n = 4. 
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Chart Young diagram Z and Admissible nmnbered Full pivot Sequence of direction Admissible pivot 
dynamical index d Young diagrams Y structure 3 for K vectors {t^i , U2, 'us, ^^4} structure J for [B, A\ 



Structure of 
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6 _l 1 1 1_ _l 2__3_^ {2,5,8,11} {62,62,62,62} {2,4,6, 
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{2,5,8,1} {61,62,62,62} {2,4,6,1} 
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{2,5,1,8} {62,61,62,62} {2,4,1,6} 
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{3,2,6,5} {62,63,62,63} {3,2,4,6} 
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{2,5,3,1} {61,63,62,62} {2,4,3,1} 



{2,5,1,3} {63,61,62,62} {2,4,1,3} 



{2,3,5,1} {61,62,63,62} {2,3,4,1} 



{2,1,5,3} {63,62,61,62} {2,1,4,3} 



{2,3,1,5} {62,61,63,62} {2,3,1,4} 



{2,1,3,5} {62,63,61,62} {2,1,3,4} 
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Table 2: Charts 6-10 for [B, A], for the case m = 3 and n = 4. 
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Chart Young diagram Z and Admissible nmnbered Full pivot Sequence of direction Admissible pivot 
dynamical index d Young diagrams Y structure 3 for K vectors {t^i , U2, 'us, ^^4} structure J for [B, A] 



Structure of 



1111 1234 



{3,6,9,12} {63,63,63,63} 



{3,4,5,6} 



+ 




(0,0,4) 



1 
















1 


1 


1 





4 
















1 


2 


3 





{3,6,9,1} {61,63,63,63} 



{3,4,6,1} 









+ * 



(1,0,3) 



3 
















1 


2 


4 





{3,6,1,9} {63,61,63,63} 



{3,4,1,5} 



{3,6,9,2} {62,63,63,63} 



{3,4,6,2} 



+ 




(0,1,3) 











3 
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2 


4 





{3,6,2,9} {63,62,63,63} 



{3,4,2,5} 
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1 
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1 







2 


4 














1 









{3,1,6,4} {61,63,61,63} 



{3, 1,4,5} 



(2,0,2) 



1 


3 














2 


4 







{1,3,4,6} {63,61,63,61} 



{1,3,4,5} 



1 
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1 







(1,1,2) 
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1 
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.■s 







2 
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4 







{3,6,1,2} {62,61,63,63} 



{3,6,2,1} {61,62,63,63} 



{3,1,6,2} {62,63,61,63} 



{3,2,6,1} {61,63,62,63} 



{3,1,2,6} {63,62,61,63} 



{3,2,1,6} {63,61,62,63} 



{3,4,1,2} 



{3,4,2,1} 



{3,1,4,2} 



{3,2,4,1} 



{3,1,2,4} 



{3,2,1,4} 



* 


* 


+ 


* 


* 


* 


* 









+ 








+ 



















+ 








* 




* 


* 


* 


+ 


* 


* 


* 


* 


* 







+ 






* 




+ 














+ 











* 


* 


* 


* 


* 


+ 


* 


* 


* 


* 


+ 









* 




* 










+ 






* 





+ 








* 


* 


* 


* 


* 


+ 




* 


* 


* 


* 


+ 











* 










+ 






* 


+ 











* 




* 


* 


* 


+ 


* 


* 


* 


* 


+ 







* 


* 


* 


* 





+ 







* 


* 


* 











+ 


* 


* 


* 


* 


* 


+ 


* 


* 


* 


* 


* 


+ 





* 


* 


* 


* 


+ 








* 


* 


* 


* 











+ 


* 


* 


* 



Table 3: Charts 1 1-15 for [B, for the case m = 3 and n = 4. 
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